Dispersal history of SARS-CoV-2 in Galicia, Spain.

The dynamics of SARS-CoV-2 transmission are influenced by a variety of factors, including social restrictions and the emergence of distinct variants. In this study, we delve into the origins and dissemination of the Alpha, Delta, and Omicron variants of concern in Galicia, northwest Spain. For this, we leveraged genomic data collected by the EPICOVIGAL Consortium and from the GISAID database, along with mobility information from other Spanish regions and foreign countries. Our analysis indicates that initial introductions during the Alpha phase were predominantly from other Spanish regions and France. However, as the pandemic progressed, introductions from Portugal and the USA became increasingly significant. Notably, Galicia's major coastal cities emerged as critical hubs for viral transmission, highlighting their role in sustaining and spreading the virus. This research emphasizes the critical role of regional connectivity in the spread of SARS-CoV-2 and offers essential insights for enhancing public health strategies and surveillance measures.

The multispecies microbial cluster of Fusobacterium, Parvimonas, Bacteroides and Faecalibacterium as a precision biomarker for colorectal cancer diagnosis.

The incidence of colorectal cancer (CRC) has increased worldwide, and early diagnosis is crucial to reduce mortality rates. Therefore, new noninvasive biomarkers for CRC are required. Recent studies have revealed an imbalance in the oral and gut microbiomes of patients with CRC, as well as impaired gut vascular barrier function. In the present study, the microbiomes of saliva, crevicular fluid, feces, and non-neoplastic and tumor intestinal tissue samples of 93 CRC patients and 30 healthy individuals without digestive disorders (non-CRC) were analyzed by 16S rRNA metabarcoding procedures. The data revealed that Parvimonas, Fusobacterium, and Bacteroides fragilis were significantly over-represented in stool samples of CRC patients, whereas Faecalibacterium and Blautia were significantly over-abundant in the non-CRC group. Moreover, the tumor samples were enriched in well-known periodontal anaerobes, including Fusobacterium, Parvimonas, Peptostreptococcus, Porphyromonas, and Prevotella. Co-occurrence patterns of these oral microorganisms were observed in the subgingival pocket and in the tumor tissues of CRC patients, where they also correlated with other gut microbes, such as Hungatella. This study provides new evidence that oral pathobionts, normally located in subgingival pockets, can migrate to the colon and probably aggregate with aerobic bacteria, forming synergistic consortia. Furthermore, we suggest that the group composed of Fusobacterium, Parvimonas, Bacteroides, and Faecalibacterium could be used to design an excellent noninvasive fecal test for the early diagnosis of CRC. The combination of these four genera would significantly improve the reliability of a discriminatory test with respect to others that use a single species as a unique CRC biomarker.

Parvimonas micra can translocate from the subgingival sulcus of the human oral cavity to colorectal adenocarcinoma.

Oral and intestinal samples from a cohort of 93 colorectal cancer (CRC) patients and 30 healthy controls (non-CRC) were collected for microbiome analysis. Saliva (28 non-CRC and 94 CRC), feces (30 non-CRC and 97 CRC), subgingival fluid (20 CRC), and tumor tissue samples (20 CRC) were used for 16S metabarcoding and/or RNA sequencing (RNAseq) approaches. A differential analysis of the abundance, performed with the ANCOM-BC package, adjusting the P-values by the Holm-Bonferroni method, revealed that Parvimonas was significantly over-represented in feces from CRC patients (P-value < 0.001) compared to healthy controls. A total of 11 Parvimonas micra isolates were obtained from the oral cavity and adenocarcinoma of CRC patients. Genome analysis identified a pair of isolates from the same patient that shared 99.2% identity, demonstrating that P. micra can translocate from the subgingival cavity to the gut. The data suggest that P. micra could migrate in a synergistic consortium with other periodontal bacteria. Metatranscriptomics confirmed that oral bacteria were more active in tumor than in non-neoplastic tissues. We suggest that P. micra could be considered as a CRC biomarker detected in non-invasive samples such as feces.

Wastewater early warning system for SARS-CoV-2 outbreaks and variants in a Coruña, Spain.

Wastewater-based epidemiology has been widely used as a cost-effective method for tracking the COVID-19 pandemic at the community level. Here we describe COVIDBENS, a wastewater surveillance program running from June 2020 to March 2022 in the wastewater treatment plant of Bens in A Coruña (Spain). The main goal of this work was to provide an effective early warning tool based in wastewater epidemiology to help in decision-making at both the social and public health levels. RT-qPCR procedures and Illumina sequencing were used to weekly monitor the viral load and to detect SARS-CoV-2 mutations in wastewater, respectively. In addition, own statistical models were applied to estimate the real number of infected people and the frequency of each emerging variant circulating in the community, which considerable improved the surveillance strategy. Our analysis detected 6 viral load waves in A Coruña with concentrations between 103 and 106 SARS-CoV-2 RNA copies/L. Our system was able to anticipate community outbreaks during the pandemic with 8-36 days in advance with respect to clinical reports and, to detect the emergence of new SARS-CoV-2 variants in A Coruña such as Alpha (B.1.1.7), Delta (B.1.617.2), and Omicron (B.1.1.529 and BA.2) in wastewater with 42, 30, and 27 days, respectively, before the health system did. Data generated here helped local authorities and health managers to give a faster and more efficient response to the pandemic situation, and also allowed important industrial companies to adapt their production to each situation. The wastewater-based epidemiology program developed in our metropolitan area of A Coruña (Spain) during the SARS-CoV-2 pandemic served as a powerful early warning system combining statistical models with mutations and viral load monitoring in wastewater over time.

ViQUF: De Novo Viral Quasispecies Reconstruction Using Unitig-Based Flow Networks.

During viral infection, intrahost mutation and recombination can lead to significant evolution, resulting in a population of viruses that harbor multiple haplotypes. The task of reconstructing these haplotypes from short-read sequencing data is called viral quasispecies assembly, and it can be categorized as a multiassembly problem. We consider the de novo version of the problem, where no reference is available. We present ViQUF, a de novo viral quasispecies assembler that addresses haplotype assembly and quantification. ViQUF obtains a first draft of the assembly graph from a de Bruijn graph. Then, solving a min-cost flow over a flow network built for each pair of adjacent vertices based on their paired-end information creates an approximate paired assembly graph with suggested frequency values as edge labels, which is the first frequency estimation. Then, original haplotypes are obtained through a greedy path reconstruction guided by a min-cost flow solution in the approximate paired assembly graph. ViQUF outputs the contigs with their frequency estimations. Results on real and simulated data show that ViQUF is at least four times faster using at most half of the memory than previous methods, while maintaining, and in some cases outperforming, the high quality of assembly and frequency estimation of overlap graph-based methodologies, which are known to be more accurate but slower than the de Bruijn graph-based approaches.

Memory-Efficient Assembly Using Flye.

In the past decade, next-generation sequencing (NGS) enabled the generation of genomic data in a cost-effective, high-throughput manner. The most recent third-generation sequencing technologies produce longer reads; however, their error rates are much higher, which complicates the assembly process. This generates time- and space- demanding long-read assemblers. Moreover, the advances in these technologies have allowed portable and real-time DNA sequencing, enabling in-field analysis. In these scenarios, it becomes crucial to have more efficient solutions that can be executed in computers or mobile devices with minimum hardware requirements. We re-implemented an existing assembler devoted for long reads, more concretely Flye, using compressed data structures. We then compare our version with the original software using real datasets, and evaluate their performance in terms of memory requirements, execution speed, and energy consumption. The assembly results are not affected, as the core of the algorithm is maintained, but the usage of advanced compact data structures leads to improvements in memory consumption that range from 22% to 47% less space, and in the processing time, which range from being on a par up to decreases of 25%. These improvements also cause reductions in energy consumption of around 3-8%, with some datasets obtaining decreases up to 26%.

Modeling the number of people infected with SARS-COV-2 from wastewater viral load in Northwest Spain.

The quantification of the SARS-CoV-2 RNA load in wastewater has emerged as a useful tool to monitor COVID-19 outbreaks in the community. This approach was implemented in the metropolitan area of A Coruña (NW Spain), where wastewater from a treatment plant was analyzed to track the epidemic dynamics in a population of 369,098 inhabitants. Viral load detected in the wastewater and the epidemiological data from A Coruña health system served as main sources for statistical models developing. Regression models described here allowed us to estimate the number of infected people (R2 = 0.9), including symptomatic and asymptomatic individuals. These models have helped to understand the real magnitude of the epidemic in a population at any given time and have been used as an effective early warning tool for predicting outbreaks in A Coruña municipality. The methodology of the present work could be used to develop a similar wastewater-based epidemiological model to track the evolution of the COVID-19 epidemic anywhere in the world where centralized water-based sanitation systems exist.

Inference of viral quasispecies with a paired de Bruijn graph.

MOTIVATION: RNA viruses exhibit a high mutation rate and thus they exist in infected cells as a population of closely related strains called viral quasispecies. The viral quasispecies assembly problem asks to characterize the quasispecies present in a sample from high-throughput sequencing data. We study the de novo version of the problem, where reference sequences of the quasispecies are not available. Current methods for assembling viral quasispecies are either based on overlap graphs or on de Bruijn graphs. Overlap graph-based methods tend to be accurate but slow, whereas de Bruijn graph-based methods are fast but less accurate. RESULTS: We present viaDBG, which is a fast and accurate de Bruijn graph-based tool for de novo assembly of viral quasispecies. We first iteratively correct sequencing errors in the reads, which allows us to use large k-mers in the de Bruijn graph. To incorporate the paired-end information in the graph, we also adapt the paired de Bruijn graph for viral quasispecies assembly. These features enable the use of long-range information in contig construction without compromising the speed of de Bruijn graph-based approaches. Our experimental results show that viaDBG is both accurate and fast, whereas previous methods are either fast or accurate but not both. In particular, viaDBG has comparable or better accuracy than SAVAGE, while being at least nine times faster. Furthermore, the speed of viaDBG is comparable to PEHaplo but viaDBG is able to retrieve also low abundance quasispecies, which are often missed by PEHaplo. AVAILABILITY AND IMPLEMENTATION: viaDBG is implemented in C++ and it is publicly available at https://bitbucket.org/bfreirec1/viadbg. All datasets used in this article are publicly available at https://bitbucket.org/bfreirec1/data-viadbg/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

From Coarse to Fine-Grained Parcellation of the Cortical Surface Using a Fiber-Bundle Atlas.

In this article, we present a hybrid method to create fine-grained parcellations of the cortical surface, from a coarse-grained parcellation according to an anatomical atlas, based on cortico-cortical connectivity. The connectivity information is obtained from segmented superficial and deep white matter bundles, according to bundle atlases, instead of the whole tractography. Thus, a direct matching between the fiber bundles and the cortical regions is obtained, avoiding the problem of finding the correspondence of the cortical parcels among subjects. Generating parcels from segmented fiber bundles can provide a good representation of the human brain connectome since they are based on bundle atlases that contain the most reproducible short and long connections found on a population of subjects. The method first processes the tractography of each subject and extracts the bundles of the atlas, based on a segmentation algorithm. Next, the intersection between the fiber bundles and the cortical mesh is calculated, to define the initial and final intersection points of each fiber. A fiber filtering is then applied to eliminate misclassified fibers, based on the anatomical definition of each bundle and the labels of Desikan-Killiany anatomical parcellation. A parcellation algorithm is then performed to create a subdivision of the anatomical regions of the cortex, which is reproducible across subjects. This step resolves the overlapping of the fiber bundle extremities over the cortical mesh within each anatomical region. For the analysis, the density of the connections and the degree of overlapping, is considered and represented with a graph. One of our parcellations, an atlas composed of 160 parcels, achieves a reproducibility across subjects of ≈0.74, based on the average Dice's coefficient between subject's connectivity matrices, rather than ≈0.73 obtained for a macro anatomical parcellation of 150 parcels. Moreover, we compared two of our parcellations with state-of-the-art atlases, finding a degree of similarity with dMRI, functional, anatomical, and multi-modal atlases. The higher similarity was found for our parcellation composed of 185 sub-parcels with another parcellation based on dMRI data from the same database, but created with a different approach, leading to 130 parcels in common based on a Dice's coefficient ≥0.5.

GeoSP: A parallel method for a cortical surface parcellation based on geodesic distance.

We present GeoSP, a parallel method that creates a parcellation of the cortical mesh based on a geodesic distance, in order to consider gyri and sulci topology. The method represents the mesh with a graph and performs a K-means clustering in parallel. It has two modes of use, by default, it performs the geodesic cortical parcellation based on the boundaries of the anatomical parcels provided by the Desikan-Killiany atlas. The other mode performs the complete parcellation of the cortex. Results for both modes and with different values for the total number of sub-parcels show homogeneous sub-parcels. Furthermore, the execution time is 82s for the whole cortex mode and 18s for the Desikan-Killiany atlas subdivision, for a parcellation into 350 sub-parcels. The proposed method will be available to the community to perform the evaluation of data-driven cortical parcellations. As an example, we compared GeoSP parcellation with Desikan-Killiany and Destrieux atlases in 50 subjects, obtaining more homogeneous parcels for GeoSP and minor differences in structural connectivity reproducibility across subjects.

Automatic group-wise whole-brain short association fiber bundle labeling based on clustering and cortical surface information.

BACKGROUND: Diffusion MRI is the preferred non-invasive in vivo modality for the study of brain white matter connections. Tractography datasets contain 3D streamlines that can be analyzed to study the main brain white matter tracts. Fiber clustering methods have been used to automatically group similar fibers into clusters. However, due to inter-subject variability and artifacts, the resulting clusters are difficult to process for finding common connections across subjects, specially for superficial white matter. METHODS: We present an automatic method for labeling of short association bundles on a group of subjects. The method is based on an intra-subject fiber clustering that generates compact fiber clusters. Posteriorly, the clusters are labeled based on the cortical connectivity of the fibers, taking as reference the Desikan-Killiany atlas, and named according to their relative position along one axis. Finally, two different strategies were applied and compared for the labeling of inter-subject bundles: a matching with the Hungarian algorithm, and a well-known fiber clustering algorithm, called QuickBundles. RESULTS: Individual labeling was executed over four subjects, with an execution time of 3.6 min. An inspection of individual labeling based on a distance measure showed good correspondence among the four tested subjects. Two inter-subject labeling were successfully implemented and applied to 20 subjects and compared using a set of distance thresholds, ranging from a conservative value of 10 mm to a moderate value of 21 mm. Hungarian algorithm led to a high correspondence, but low reproducibility for all the thresholds, with 96 s of execution time. QuickBundles led to better correspondence, reproducibility and short execution time of 9 s. Hence, the whole processing for the inter-subject labeling over 20 subjects takes 1.17 h. CONCLUSION: We implemented a method for the automatic labeling of short bundles in individuals, based on an intra-subject clustering and the connectivity of the clusters with the cortex. The labels provide useful information for the visualization and analysis of individual connections, which is very difficult without any additional information. Furthermore, we provide two fast inter-subject bundle labeling methods. The obtained clusters could be used for performing manual or automatic connectivity analysis in individuals or across subjects.

Efficient processing of raster and vector data.

In this work, we propose a framework to store and manage spatial data, which includes new efficient algorithms to perform operations accepting as input a raster dataset and a vector dataset. More concretely, we present algorithms for solving a spatial join between a raster and a vector dataset imposing a restriction on the values of the cells of the raster; and an algorithm for retrieving K objects of a vector dataset that overlap cells of a raster dataset, such that the K objects are those overlapping the highest (or lowest) cell values among all objects. The raster data is stored using a compact data structure, which can directly manipulate compressed data without the need for prior decompression. This leads to better running times and lower memory consumption. In our experimental evaluation comparing our solution to other baselines, we obtain the best space/time trade-offs.